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ABSTRACT 

We present the results of a gravitational lensing analysis of the bright = 2.957 sub- millimeter 
galaxy (SMG), HERMES J105751. 1+573027 found in the iJer^c/ie// SPIRE Science Demonstration 
Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high res- 
olution imaging available in optical and Near-IR channels, along with CO emission obtained with 
the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and 
hence estimate the total lensing magnification to be /i = 10.9 ± 0.7. We measure the half-light radius 
i^eff of the source in the rest-frame Near-UV and V bands that characterize the unobscured light 
coming from stars and find i^eff,* = [2.0 ±0.1] kpc, in good agreement with recent studies on the Sub- 
Millimeter Galaxies population. This lens model is also used to estimate the size of the gas distribution 
(i?eff,gas = [1.1 ± 0.5] kpc) by mapping back in the source plane the CO ( J = 5 ^ 4) transition line 
emission. The lens modeling yields a relatively large Einstein radius i^Ein = 4V10±0V02, corresponding 
to a defiector velocity dispersion of [483±16] kms"^ This shows that HERMES J105751. 1+573027 is 
lensed by a galaxy group- size dark matter halo at redshift z\ ~ 0.6. The projected dark matter contri- 
bution largely dominates the mass budget within the Einstein radius with /dm(< ^Ein) ^ 80%. This 
fraction reduces to /dm(< ^eff,Gi — 4.5 kpc) ~ 47% within the effective radius of the main defiecting 
galaxy of stellar mass M^^gi = [8-5 ±1.6] x lO^^M©. At this smaller scale the dark matter fraction is 
consistent with results already found for massive lensing ellipticals at z 0.2 from the SLAGS survey. 
Subject headings: gravitational lensing - submillimeter - galaxies: elliptical and lenticular, cD - 
galaxies: halos 
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1. INTRODUCTION 



The current generation of wide field surveys at sub- 
millimeter and millimeter wavelengths is now providing 
us with large numbers of high redshift galaxies containing 
large amounts of dust heated by intense star formation 
or AGN activity (e-g-, Hughes et al.||1998| iBarger et al. 
1998; Blain et al.||2QQ2rTChapman et al.|[2QQ5; Coppin 



et al.||2QQ6 Austermann et a l.||2QlQ ). This population 
ot sub-millimeter galaxies (SMGs) is easily detectable in 
the redshift range 1 ^ ^ ^ 5 thanks to a strong negative 
k-correction when observed at wavelengths A > 500 jam. 
This property, along with a sharp fall-off at the bright lu- 
minosity end of their luminosity function, makes bright 
SMGs good candidates for being strongly gravitation- 
ally lensed fe.g., | Blain 1996; Negre llo et al. 2007; CoorayJ 
|et al. 2010!lNegrello et al. 2010; V ieira et al., 2010J .~Er 
ncient identification of lensed SMGs is only now becom- 
ing possible thanks to surveys like the Herschel Multi- 
tiered Extragalactic Survey (HerMES) or the Herschel - 
Astrophysical Terahertz Large Area Surve yfH- ATLAS) 



conducted with the Herschel satellite (e 



2010 Griffin et al. 2010 Negrello et al. 20 10D and tfom 



the g round with the South Pole 'I'elescope ( | Vieira et al 



Oliver et al 

I 

Vi 



2010|). 

The interest of building large samples of lensed SMGs 
for getting better insights on the properties of these oth- 
erwise very faint objects is clear, and recent results are 
already shedding some light on the spatial d istribution 



of gas, dust and stars in these SMGs (e.g., Swinbank 
et al.|201Qb ). However the redshift distribution of SMGs. 



which peaks in the ra nge 2 — 2.5 ( Chapman et al.||2005 
[Wardlow et al. 2010[ ), is also well suited to probe the 
mass distribution of high 



redshift deflectors, typically in 
the range 0.3 < zi < 1.5, which complements local stud- 
ies like the Sloan Lens ACS Survey (SLAGS) which ar e 
limited to z < 0.4 ( [Bolton et al.|2008| [Auger et al.|2010"l ). 

In this paper, we present the modeling of the gravita- 
tionally lensed SMC HERMES J105751. 1+573027, also 
referred to as HLSW-01, discovered in Herschel/ SFIRE 
observations during the Science Demonstration Phase 
by the HerMES program. Its J2000 coordinates are 
RA=10:57:51.0, Dec=+57:30:25 .8 with a lens photomet- 
ric redshift of zi = 0.60 ± 0.04 (|Oyaizu et al.||2008) and 
a source redshift of z^ = 2.9575 ± 0.0001 (Riechers et al. 
2011; Scott et al. 2011; hereafter Rll and Sll). The 
goal of this lens modeling work is two-fold. First we 
want to recover the intrinsic light distribution of the 
source, while optimally taking advantage of the magnify- 
ing power of the deflector. Detailed investigations on the 
lensed source are developed in associated papers (Conley 
et al. 2011 [hereafter Gil], Rll, Sll). And second we 
want to probe the mass content of the foreground object 
which, given the large image separation of the multiple 
images, might be very massive. 

The paper is thus organized as follows. In ^we present 
the lens modeling techniques and the optical and Near- 
IR data we shall use along with the main results. The 
best fit lens model is then used in ^to reconstruct the 
source GO (J = 5^4) light distribution at 576 GHz ob- 
served with the Plateau de Bure Interferometer (PdBI). 
In Q we interpret the lens model results to measure the 
balance of dark and luminous matter in the inner 30 kpc 
of the deflector. We conclude in ^ Throughout, we 



assume a concordance cosmology with matter and dark 
energy density ft^ = 0.3, I^a = 0.7, and Hubble constant 
Hq=70 km s~^Mpc~^. Magnitudes are expressed in the 
AB system. 

2. LENS MODELING OF OPTICAL DATA 
2.1. Observations 

For accurate lens modeling, we use the best spatial 
resolution images of sufficient signal-to-noise that are 
currently available for HLSW-01. A 1400 second image 
of the system was taken using Laser Guide Star Adap- 
tive Optics in the band with the NIRG2 instrument 
mounted on the Keck II telescop^ The observing con- 
ditions allowed us to achieve a typical 0^2 FWHM point 
spread function. The FWHM is well sampled with a 
plate scale of 0V02. We additionally use a 1 hour Subaru 
SuprimeGam i band observation with ~ 0^.^74 FWHM 
seeing and 0^^202 pixel size. 

A Kp image of the central region of the system is shown 
in Fig. [T] The lensing configuration can easily be de- 
scribed despite the relatively low surface brightness of 
the multiply-imaged features of HLSW-01. The config- 
uration of multiple images lifts any ambiguities as we 
clearly see four images of similar surface brightness form- 
ing a, so-called, fold configuration with images 1 and 2 
presumably merging through the critical line and images 
3 an d 4 being lower magni fication conjugate images (see 
e.g., [Schneider et "aLl|1992[ for a description of catastro- 
phe theory in the context of gravitational lensing). A 
closer look at the band image shows that image 1 is 
perturbed by a small galaxy, G4, which seems to be mas- 
sive enough to split image 1 into two pieces (la and lb) 
on both sides of G4. Galaxies G2 and G3 may also act 
as potential perturbers on the innermost multiple image, 
4. To a lesser extent, G5 might also be considered as a 
perturbing galaxy. Gl is the central galaxy of the mas- 
sive deflecting structure, which presumably should be a 
group of galaxies, given the large S" image sepa- 

ration. Having no redshift information at all for these 
perturbers, the simplest assumption is to consider them 
to be at the same redshift as Gl, keeping in mind that 
changes from this hypothesis are of second order in the 
lens modeling. Indeed, the total magnification of HLSW- 
01 is not changed, while the absolute mass calibration of 
G2...G5 depends on the assumption that they are at the 
same redshift as Gl. However, measuring the mass of 
perturbers is not the main motivation of the modeling. 

2.2. Method 

The lens modeling builds on the dedicated code sl_f it 
previously used fo r gal axy-scale stro ng lenses (e.g., 
Gavazzi et al ||2007| [2QQ8| [Ruff et al.||201 ) . It fits model 
parameters oi simple analytic lensing potentials. sl_fit 
can be run in three different regimes of increasing com- 
putational cost. The first mode makes use of the coordi- 
nates of image plane points and minimizes the distance 
to their parent source plane locations in a way similar 
to gravlen s ( Keeton 2001) or lenstool (Kneib 1993 
Julio et al.J[2007|). The second mode uses the lull surface 



brightness distribution and attempts to account for it 

^ |http : //www2 . keck . hawaii . edu/realpublic/inst/nirc2/| 

I index .html I 
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Fig. 1. — Xp band adaptive optics-corrected overview of HLSW- 
01 with labeled multiple images (red) and foreground galaxies 
(black). North is up and East is left. The bottom right inset zooms 
into the image 1 that appears to be perturbed by the foreground 
galaxy G4 and split into two pieces. 

with one or more simple analytic light components that 
we take to have a unique Gaussian radial profile with el- 
liptical shape (see e.g., Marshall et al.|2007 Bolton et al. 
|2QQ8[ for similar techniques), finally the tmrd mode im- 
plements a pixellized linear reconstruction of the source 
plane light distributi on while fitting for the non-linear 



potential parameters fWarren fc Dye|2QQ3 Treu fc Koop 
[mans 2004; Suyu et al. 2006); we did not consider the 



latter mode here as its computational cost is prohibitive 
for the large images and complex gravitational potential 
of this system. 

The lensing potential is assumed to be made of a 
cored isothermal ellipsoid, centered on the main deflector 
galaxy Gl, and that is supposed to capture the lensing 
contribution of the dark matter halo as well as the stel- 
lar component of Gl. Given the absence of a radial arc 
or ce ntral demagnified images (see e.g., 'Gavazzi et al. 



2003[ ), the details of the assumed potential in the inner- 
most parts (r < 2'^) of the lens should not be important. 
The peak of Gl's light distribution is assumed to be the 
center of this potential component. The convergence pro- 
file of the central mass component is given by 



t(0 = 



1 



Ve 



(1) 



where the scaling parameter 6cent is related to the 
velocity dispersion of the deflector through b/V^ = 
(cr^/186.21 kms"^)2As/A, the core radius is Tc, ^ = 
jc^ is the radial coordinate that accounts for the 
ellipsoidal symmetry of the isodensity contours and q is 
the minor-to-major axis ratio. The orientation of the ma- 
jor axis, ^cent7 relative to the x-axis is allowed to vary, 
although this is not explicit in the definition of ^. As 
will be seen below, lens modeling of extended images 
can yield formal errors on h of order one percent, and 
therefore similar errors on cr^. However, here we propa- 



gate uncertainties in D\^/D^ due to the relatively poorly 
known lens redshift {z\ = 0.60 ± 0.04). This results in 
a dominant additional 3% error that we add in quadra- 
ture to model uncertainties on 6, keeping in mind that 
a spectroscopic measurement of zi would readily reduce 
lens modeling errors on ay to the percent level. 

We carefully account for the perturbing galaxies G2, 
G3, G4 and the more distant G5, and model them as 
truncated isothermal ellipsoids, each with a free velocity 
dispersion and a constant cutoff radius that we fix to 
a value of rt = 1^5 corresponding to lOkpc (see e.g., 
[Limousin Tt al.||2007H2009| |Suyu fc Halkola||2010|). For 



each ot these galaxies the convergence protile takes the 
form: 



^pert 



{r- rG,i) 



^pert,i 



1 



1 



2 



(2) 



SO that the total mass of a perturber is Mpert,i = 
7rScrit^pert,i^t whcrc Scrit IS the Critical surface density, 
which for this system has a value Scrit — [1.80 ± 0.05] x 
10^ Mo kpc"^ [8.05 ±0.24] x 10^° M© arcsec-^. The 
3% quoted errors are dominated by the uncertainty in the 
photometric redshift of the lens, but this is still negligible 
compared to the uncertainty in the velocity dispersion of 
the main component. As an approximation, we match 
the ellipticity and orientation of these perturbing sys- 
tems to that of the stars they host, although we note 
that recent analyses have shown that departures from 
this simple assumption might occur and is sometimes 
observed with more favorable lensing config^urat ions and 
deep space-based optical datafSuyu fc Halkola||2010 ). 

To summarize, the lensing potential is described by 
4 parameters for the central component: 6cent, ^c, 
^cent, Ocent and wc also allow the 4 scaling parameters 
6pert,i=2...5 of the perturbing potentials to vary. 

2.3. Preliminary analysis: conjugate points 

We first attempted to fit the model parameters of 
the main central component assuming that the source 
is point-like (first mode of sl_fit) and the perturbers 
have a negligible effect. By ignoring the effect of G4 we 
had to neglect the additionally produced image lb and 
just considered the brightest pixel of each of the images 
la, 2, 3 and 4 as input constraints on the potential. We 
performed a fit to the parameters Qcent^ ^cent, ^cent and Tc- 
We used these results to further understand the source 
structure (fainter spots) seen in image 2 and how it is cast 
onto other images la, 3 and 4 respectively, and, hence in- 
crease the number of constraints (see e.g., Gavazzi et al 
|2003[ for an illustration of the process). I'his also al- 
lowed us to determine which part of the source should 
be inside the extra caustic caused by the perturber G4 
and thus seen as image lb. Altogether we identified 4 
multiply-imaged knots. 

We ran a Monte-Carlo Markov Chain sampler of the 
posterior distribution of the 8 model parameters related 
to the gravitational potential. At this stage the most 
important result^ are that we find quite an elongated 
mass distribution with an axis ratio ^cent = 0.41 ± 0.03 

^ For a given parameter we quote the median and 68% confi- 
dence level errors as given by the 16th and 84th percentiles of the 
marginalized distribution. 
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Fig. 2. — Results of the surface brightness fit of the Keck LGSAO 
Kp band imaging. Top left panel: Input image with the foreground 
defiectors subtracted off (with their core completely masked out in 
some cases) using galf it. Top right panel: Image plane model pre- 
dictions. Bottom left panel: Image plane residuals (data-model). 
Bottom right panel: Source plane model predictions. In this latter 
panel, the caustic lines are shown in red, whereas the other panels 
show the image plane critical lines. Scales are given in arc seconds. 

and orientation ^cent = 13?8 ± 0?5 (North to East coun- 
terclockwise). The circularized Einstein radius is i^Ein = 
4V02 ± 0V05, corresponding to a velocity dispersion of 
ay = [473 ± 15] kms~"^, which suggests that the deflec- 
tor is a massive group of galaxies. The modeling seems 
to require a finite core radius Tc = [2.0 ± 1.0] kpc, but 
it is difficult to interpret this as a constraint on the 
dark matter distribution since the central component 
of Eq. ([T]) captures both the contribution of stars and 
dark matter. The mass of perturber G4 is relatively 
weh constrained: Mqa = [25 ± 5] x lO^^M©, whereas 
Mg2 = [13 ±9] X lO^^Mo, Mg3 = [64 ±6] x lO^^M© and 
Mg5 = [20 ± 12] X lO^^M©. G5 has a nearly negligible 
effect on the lensing configuration. 

This best fit model predicts local magnification^ of 
-1.87, 4.37, 2.57 and -1.86 for images lb, 2, 3 and 4, 
respectively, which means that the total magnification 
should be of order 10 for a point-like source. However 
the above treatment of point-like lensing observables, al- 
though it allowed us to quickly explore the space of lens 
potential parameters, does not take full advantage of the 
spatially resolved surface brightness distribution. It thus 
cannot give us a clear idea of the intrinsic source light 
distribution, on which the actual magnification factor de- 
pends. 

2.4. Keck LGSAO hand light distribution 

We therefore placed an extended source with an an- 
alytic light profile made of an elliptical Gaussian with 
free central position (xg, ^g), axis ratio (g'g), position an- 
gle 6>g, fiux Fg and half-light radius i^efT,s- We optimized 

^ Negative values indicate a change of image parity with respect 
to that of the source. 
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Fig. 3. — Change of the total magnification as a function of the 
characteristic source size i^efF,s- The error bar shows the statistical 
error inferred from fitting the Kp band Keck LGSAO image (cor- 
responding to a rest frame V band observation). The drop at low 
values of i?eff,s is artificially due to sampling limitations but does 
not affect the analysis. 

these parameters along with the potential parameters al- 
ready introduced in the previous exploration phase us- 
ing the surface-brightness mode of sl_fit. Fig. [2] shows 
the resulting fit in the image plane. The correspond- 
ing reconstructed source is shown in the bottom right 
panel. The best fit is achieved with an extended source 
of i^p band magnitude 22.12 ± 0.08 and half-light radius 
^eff,s = [1.88 ±0.09] kpc. 

We estimate the net magnification experienced by the 
source through a direct numerical sum of pixel values in 
the image and the source planes and find the magnifica- 
tion to be /i = 10.86±0.68 after marginalizing over source 
and potential parameters. We show in Fig. |3] the change 
of the magnification as a function of the source size i^eff,s- 
This suggests that, provided the size stays within a factor 
of ~ 2 for one band to another, the systematic change in 
magnification remains within the current statistical er- 
rors. We can also see that magnification is shown to be 
< 11.5 whatever the source size and this can readily be 
cast as a robust lower bound on the intrinsic source fiux 
for a given total observed flux. 

Compared to the previous analysis using the conjuga- 
tion of bright knots, the fit of an extended source yields 
little change in the recovered lens potential. The key 
features of the lens potential are 

• The axis ratio of the mass distribution is ^cent = 
0.45 ± 0.02, with orientation l9cent = 11°9 ± 0?4. 

• The (circularized) Einstein radius is i^Ein = 4^10 =b 
0^.^02, which corresponds to a velocity dispersion 
of (jy = [483 ± 16] kms"^, in agreement with the 
preliminary study of §2.3[ The mass content of this 
defiecting component is further discussed in Q 

• The core radius Tc = [2.7 ± 0.8] kpc is found to be 
slightly larger. It is mostly constrained by the size 
and shape of the innermost image 4. 

• There is little change in the mass of perturbers, 
for which we find: Mqa = [26.8 ± 1.8] x lO^^M©, 
Mg2 = [3.2 ± 2.4] X lO^^Mo, Mgs = [49.4 ± 5.2] x 
lO^^M© and Mg5 = [12.7 ± 7.0] x lO^^M©. 
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Fig. 4. — Subaru i band modeling results. Top left panel: input 
image with the foreground deflectors subtracted off. Top right 
panel: image plane model predictions. Bottom left panel: Image 
plane residuals (data-model). Bottom right panel: source plane 
model predictions. 

2.5. Seeing -limited Subaru i band light distribution 

We now take advantage of the somewhat deeper seeing 
hmited Subaru i band image to investigate the abihty 
of the Kp model to account for observations at shorter 
wavelengths. This image corresponds to rest frame ~ 
2000A NUV emission that we expect to be lumpier and 
severely obscured by dust, and thus not to have the same 
extent as the rest frame V band probed by the Kp band 
image. 

We therefore repeated the previous analysis on the 
Subaru i image, but we considered the best fit potential 
values above and only attempted to fit for the i band 
light profile parameters of the source: central position 
(xg, ^s)^ axis ratio (qs)^ position angle 0^^ flux and 
half-light radius i^eff,s- Fig.[4]shows the result of the sur- 
face brightness distribution nt in the image plane. The 
corresponding reconstructed source is shown in the bot- 
tom right panel. The source is found to have a similar 
appearance in i and bands with a half-light radius 
^eff,s = [2.00 ± 0.01] kpc and an i band AB magnitude 
of 22.72 i 0.01. These small statistical errors should be 
treated with caution, as but the residuals are worse 
than the for case, which might indicate a more com- 
plex intrinsic NUV light distribution. These departures 
from a simple Gaussian elliptical profile could not be 
observed in because of the substantially lower signal- 
to- noise ratio. 

This latter analysis demonstrates that lens model pa- 
rameters found by fitting the light distribution in the 
redder Kp filter give satisfying results for the fit of the 
light distribution in a different band, i. 

3. CO(5-4) LINE DISTRIBUTION 

Given the above success at accounting for the lensed 
visible/NIR light distribution of HLSW-01, we extend 
our analysis to sub-mm wavelengths using spatially re- 



solved PdBI observations of the CO (J = 5 4) transi- 
tion line at 576 GHz (see Rll for details). 

Since our lens modeling code was originally designed 
to model optical/NIR images, we could not fit these ob- 
servations directly in the uv-plane, which would allow 
a more careful handling of the beam convolution and 
limit the effect of noise correlations. In order to circum- 
vent these difficulties we used a reduced image that was 
CLEANed with a synthesized beam of 4^8 x 2^7 FWHM, 
with a major axis oriented +51? 6 East of North. The 
measured noise rms is 1.2 mJy beam" ^. 

We did not attempt to fit for lens potential model 
parameters, as they were better determined with visi- 
ble/NIR data, but we fit for source position integrated 
intrinsic flux and shape. Fig. |5] shows the result of the 
CO (J = 5 4) emission map fit in the image plane. 
The corresponding reconstructed source is shown in the 
bottom right panel. We can see that the modeling yields 
very small residuals beyond the secondary lobes that can- 
not be captured by our direct space modeling strategy. 
Our main inferences for the source parameters are 

• We find the source half flux radius to be i^eff,s = 
[1.13 ± 0.53] kpc, slightly smaller than the extent 
we found for the rest-frame V and NUV bands fits. 

• We measure a small offset between the center of the 
gas distribution as probed by the CO (J = 5^4) 
transition and the stars that emit at visible/NIR 
wavelengths. The typical offset is ARA = — 0^.^29 ± 
0'.'07 and ADec = -OVlO ± 0V05 (corresponding to 
about 2.4 kpc). This is suggestive that rest frame 
UV and visible light is severely obscured by the 
large dust content presumably associated with the 
gas revealed by the CO(J = 5^4) line (Rll). 
The presence of dust is confirmed by the Far IR 
emission (Cll). We stress however that further in- 
vestigation of this is left for future work, including 
higher resolution mm or radio data. 

• The intrinsic source luminosity of the CO (J = 5 

4) line is i^co(5-4) = [3.76 ±0.44] x 10^° K km s"^ 
pc^ including the equally important instrumental 
and modeling error terms. 

The good signal-to-noise ratio of the PdBI data sug- 
gests that one could use the kinematical information to 
try and reconstruct the intrinsic velocity field as traced 
by the CO (J = 5^4) transition line. However, because 
of insufficient spatial resolution and perhaps non trivial 
patterns in the velocity field, we were not able to explain 
the marginally resolved velocity field presented in Fig. 5 
of Rll. 

4. MASS CONTENT OF THE DEFLECTOR 

We turn now to an analysis of our findings regarding 
the lens potential and compare the inferred mass distri- 
bution with the propert ies of the main deflector. 

We used galfit (Peng et al.'"2002) to measure the 
Kp light properties of the main galaxy Gl and found 
that th e surface bri ghtness profile is well fit by a Sersic 
profile ( [Sersic 1968) of index n = 3.06 ± 0.10, effective 
radius i?eff,Gi = [4.5 ± 0.1] kpc, total magnitud^ 

^ Photometry is quite uncertain with adaptive optics imaging, 
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stantially different from lensing findings i n z ^ 0.2 el- 





FlG. 5. — Plateau de Bure Interferometer CO(J = 5—^4) line 
flux density maps modeling results. Top left panel: input image 
with the foreground deflectors subtracted off. Top right panel: 
image plane model predictions. Bottom left panel: Image plane 
residuals (data-model). Bottom right panel: source plane model 
predictions. Note the slightly larger scale compared to Figs. |2] and 

17.41 ± 0.10 and a magnitude 18.37 ± 0.04 in the i band. 
Assuming a spectral energy distribution typical of an el- 
liptical galaxy (Coleman et al.||198 0), we can transform 
this observed magnitude mto a rest-frame V band lu- 
minosity Lv = [3.9 ± 0.6] X 10^^ L0, including errors on 
photometry (4%), photometric redshift (16%) and uncer- 
tainties in filter conversions (13%). Taking into account 
the luminosity decline due to t he passive evolution of its 
old stellar populations ( Treu et al. 2001) , 



d M* 

log _^ - -0.40 ± 0.05 , 

dz Ly 



(3) 



we can predict the evolution-corrected luminosity. Ap- 
ply ing the same stellar mass-to-light ratio as found 
by Gavazzi et al. (2007) for massive early- type galax- 
ies in the SLACS survey we get a total stellar mass 
M*,Gi = 8.5 ± 1.6 X lO^^M© for the central galaxy. 
This large stellar mass, along with the large charac- 
teristic velocity dispersion inferred from lens modeling 
(jy = [483 ± 16] kms""*^, suggests that the deflector is 
likely a massive group of galaxies with Gl being its dom- 
inating central galax}{^ 

It is tempting to compare these results with the SLAGS 
findings for massive Early- Type Galaxies at redshift z ~ 
0.2 and with median stellar mass ~ 2 — 4 x 10^^ (de- 
p ending on the ch oice of the stellar Initial Mass Function 
Auger et al.||2010 ). We find that our lens model requires 
a finite core radius Tc = [2.72 ± 0.84] kpc. This is sub- 

and calibration was made difficult due to uncertainties in conver- 
sion between NIRC2 Kp filter and the reference 2MASS Ks pho- 
tometry. 

^ We could not find any concentration of galaxies within several 
arcmin probed by the Subaru image, further indicating that Gl is 
the main mass component. 



liptical galaxies of slightly lo wer mass (e.g., Bolton et al. 
2008[ Koopmans et al.|2009 ) in which a Singular Isother- 
mal Elliptical profile provides a good description of the 
total lensing potential. In addition the lensing-inferred 
velocity dispersion is unlikely to reflect the stellar veloc- 
ity dispersion of Gl that would be measured by opti- 
cal/NIR spectroscopy. Indeed, the Fundamental Plane 



(FP) of early- typ e galaxies (Djorgovski & Davis |1U^Y 
Faber et al. 1198 71) predicts that the velocity dispersion o: 



1987 



Gl sh ould be [400^30] km s~^ with the latest normaliza- 
tion of Auger et al. ( |2010 ). However the finite core radius 
requires a numerical integration of the Jeans equation to 
relate our total density profile inferred from lensing and a 
true velocity dispersion dap within a fiducial aperture of 
radius i^ao = ^'-^ — 10 kpc. Following the prescriptions 



of Mamon & Lokas ( 2005| ) for fast numerical integrations 
of the Jeans equation, assuming an isotropic velocity ten- 
sor, our model predicts a value dap = [380 ± 5] kms~^, 
in good agreement with the FP value. 

A direct comparison with SLAGS is complicated by the 
very different values of Einstein radii. For SLAGS the ra- 
tio i^Ein/^eff ^ 0-5 whcrcas we observe i^Ein ^ 27 kpc ~ 
6 X i^eff- Therefore, it is not surprising that the total 
mass within i^Ein is a factor of ~ 5 greater than the stel- 
lar mas^ The deflector is thus a group of galaxies that is 
more dark matter dominated than most SLAGS massive 
early-type galaxies. We note that the cosmological crit- 
ical density for lensing Ecrit = c^/(47rG)I)s/(Ms^os) is 
more favorable for SMGs (with typical redshifts Zs ~ 2 — 3 
and deflectors at zi ~ 0.6) than for SLAGS lenses (for 
which Zs ~ 0.7 and zi ~ 0.2). Gonsequently, the de- 
flector in HLSW-01 and similar z ^ 0.6 systems will 
typically have twice as large an Einstein radius as lower 
redshift SLAGS lenses, regardless of differences in the de- 
flector mass. This large radius implies that these lenses 
will presumably be more dark matter dominated. Ex- 
trapolating our best fit density profile inward to the 
effective radius yields a projected dark matter fraction 
/dm(< i^eff) = 0.47^0.26. agreement with SLAGS find- 
ings. 

5. CONCLUSIONS 

We have calculated a detailed lens model of the newly 
found gravitational lens HERMES J105751. 1+573027. 
Taking advantage of the best available spatial resolution 
images in the i^p-band using the Laser Guide Star Adap- 
tive Optics system at the Keck II telescope, in the i-band 
with the Subaru telescope and using Plateau de Bure In- 
terferometer to probe the GO (J = 5^4) transition line 
at 576 GHz, we were able to infer the mass distribution 
in the inner ~ 30 kpc of the deflecting structure, which 
turns out to be a massive galaxy, presumably at the cen- 
ter of a galaxy group size halo of characteristic velocity 
dispersion = [483 ± 16] kms~^. 

The redshift distribution of lensed SMGs will natu- 
rally select objects in a favorable range for producing 
large Einstein radii, compared to local lenses. This pro- 
vides a great opportunity to probe the total density pro- 
file and the dark matter content of intermediate redshift 

^ This corresponds to a total mass-to-light ratio of Mtot/Ly = 
18-0^2 7 (^Z^^)©' cigain, correcting for luminosity evolution 
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halos with high accuracy. The lensing configuration of 
HLSW-01 yields a value of i?Ein = 4'/05 - 27kpc that 
corresponds to 6 times the effective radius of the cen- 
tral galaxy. At these scales we observe a projected dark 
matter fraction of about 80%. The defiector of HLSW- 
01 probably stands in an interesting transition regime 
between cluster and galaxy scale lenses, in which a joint 



lensing plus dynamical analysis ( Miralda-Escude 



Sand et al. 2004 Gavazzi 2005 Newman et al. 



would be very informative tor the actual small-scale dar 
matter distribution. 

Much of the novelty of the large number of lensed 
SMGs to be uncovered in ongoing and future submilli- 
metric surveys resides in the opportunity for studying 
faint heavily star-forming galaxies with good resolution 
and signal-to-noise that would otherwise be unreach- 
able. In the particular case of HLSW-01, the source 
turns out to be magnified by a factor ~ 10, which al- 
lowed us to measure the extent of the emitting gas in 
the CO (J = 5 4) transition line as well as young 
stars emitting in the NUV and visible bands. We mea- 
sured a half-light radius i^eff,gas = [1.13 ± 0.53] kpc and 
^eff,* = [2.0 ± 0.1] kpc using PdBI and Keck data, 
respectively, in ex cellent agreement with th e sample of 



SMGs studied by [Swinbank et al | ( |2010aD . There is 



some evidence that the peak of the CO emission and 
the peak of the NUV/ visible light could be offset by 
~ 0^3 ^ 2.4 kpc, which could be explained if the stellar 
light is obscured by dust associated with the cold gas. 
The source reconstruction allows for a careful estimate 
of the intrinsic source fiux at many other wavelengths. 



thus allowing a more detailed analysis of the gas and 
dust content of this peculiar SMC (Cll, Rll, Sll) and 
its gas content. 

The large number of lensed SMGs like HLSW-01 to 
be found in the coming years with HerMES, H- ATLAS, 
the SPT and other surveys, will allow detailed analyses 
of the central regions of massive galaxies with exquisite 
accuracy over a broad range of defiector redshifts. This 
is a good opportunity to constrain evolutionary trends in 
their dynamical properties, including clues on the role of 
dry vs. wet mergers or the role of the central AGN. At the 
same time, the resolving power of gravitational lensing 
will allow detailed investigations of the stellar, gaseous 
and dust content of massively star forming SMGs at high 
redshift. 



SPIRE has been developed by a consortium of in- 
stitutes led by Cardiff Univ. (UK) and including 
Univ. Lethbridge (Canada); NAOC (China); CEA, LAM 
(France); IFSI, Univ. Padua (Italy); lAC (Spain); Stock- 
holm Observatory (Sweden); Imperial College London, 
RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); Cal- 
tech, JPL, NHSC, Univ. Colorado (USA). This develop- 
ment has been supported by national funding agencies: 
CSA (Canada); NAOC (China); CEA, ONES, CNRS 
(France); ASI (Italy); MCINN (Spain); SNSB (Sweden); 
STFC (UK); and NASA (USA). The IRAM Plateau 
de Bure Interferometer is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 
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